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^ , Abstract 

1^ i I review the analysis of (2+l)-dimensional Yang-Mills {YM2+1) theory 

via the use of gauge-invariant matrix variables. The vacuum wavefunction, 
string tension, the propagator mass for gluons, its relation to the magnetic 
mass for YM^^i at nonzero temperature and the extension of our analysis to 
^ . the Yang-Mills-Chern-Simons theory are discussed. A possible extension to 

! 3 + 1 dimensions is also briefly considered. 

O 

5 : 1 Introduction 



In this talk I shall review a Hamiltonian approach to Yang-Mills theory in 



D ! two spatial dimensions {YM2+1), on which I have been working for the last 

^ ' few years. In this analysis, it is possible to carry out some nonperturbative 

calculations to the extent that results on mass gap and string tension can be 
^ ■ compared with lattice simulations of the theory. I shall also discuss recent 

■ attempts to extend this to 3 + 1 dimensions. The work I shall report on was 

done in collaboration with D. Karabali and Chanju Kim [1, 2, 3] and more 
recently with A. Yelnikov [4]. 

Let me begin with a few general comments. The 2 + 1 dimensional case 
is interesting because it can be a model simple enough to analyze mathe- 
matically and yet nontrivial enough to teach us some lessons about (3+1)- 
dimensional YM theories. Another reason to study YM2+1 is its relevance 
to magnetic screening in YM^^i at high temperature. Gauge theories at 
finite temperature have worse infrared problems than at zero temperature 

*Talk at "Space-time and Fundamental Interactions: Quantum Aspects", conference 
in honor of A. P. Balachandran's 65th birthday, Vietri sul Mare, Salerno, Italy, May 2003; 
to be published in Mod. Phys. Lett. 



1 



due to the divergent nature of the Bose distribution for low energy modes. 
A dynamically generated Debye-type screening mass will eliminate some of 
these, but we need a magnetic screening mass as well to have a perturbative 
expansion which is well defined in the infrared. The mass gap of the 2 + 1 
dimensional theory at zero temperature can be interpreted as the magnetic 
mass of the 3 + 1 dimensional theory at high temperatures, with a certain 
relationship between the coupling constants [5]. As for the 3 + 1 dimensional 
theory at zero temperature, being part of the standard model, its importance 
needs no emphasis. 

For YM2+1 the following observations will be useful. The coupling con- 
stant has the dimension of mass and it docs not run as the 3+1-dimensional 
coupling does. The dimensionless expansion parameter of the theory is k/e"^ 
or e^/k, where k is a. typical momentum. Thus modes of low momenta 
must be treated nonpcrturbativcly, while modes of high momenta can be 
treated perturbativcly. There is no simple dimensionless expansion param- 
eter. YM2+1 is perturbativcly super-rcnormalizable, so the ultraviolet sin- 
gularities are well under control. The 3 + 1 dimensional theory will require 
renormalization and the associated dimensional transmutation. 

2 The 2 + 1 dimensional theory 

We shall first discuss YM2+1. Consider a gaiigc theory with group G = 
SU{N) in the Aq = gauge. The gauge potential can be written as Ai = 
—it°'A°^, i = 1,2, where t"' are hermitian N x A^-matrices which form a basis 
of the Lie algebra of SU{N) with = if'H^, Tr(rt^) = ^S^K The 

spatial coordinates Xi,X2 will be combined into the complex combinations 
z = Xi — 1x2, z = Xi + 1x2 with the corresponding components for the 
potential A = A^ = \{Ai + iAg), A = A~^ = \{Ai - 2^2) = -{A^^ . The 
starting point of our analysis is a change of variables given by 

A^ = -d^MM~\ A^ = M^-^d^M^ (1) 

Here M, are complex matrices in general, not unitary. If they are unitary, 
the potential is a pure gauge. The parametrization (1) is standard in many 
discussions of two-dimensional gauge fields. A particular advantage of this 
parametrization is the way gauge transformations are realized. A gauge 
transformation A^ — >• A\^^ = g^^Aig + g'^dig, g{x) G SU{N) is obtained 
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by the transformation M M^^^ = gM. The gauge-invariant degrees of 
freedom are parametrized by the hermitian matrix H = M'^M. Physical 
state wavefunctions are functions of H. 

There will be three basic steps in our analysis, namely, the determination 
of the inner product of the wavefunctions by evaluating the gauge-invariant 
volume measure, rewriting the Hamiltonian as an operator involving the new 
variables and finally using these to solve the Schrodinger equation. We will 
consider these in turn. 

The measure of integration over the fields A, A is dii{A)/vol{Q^) where 
d^{A) = Yl,j. adA°'{x)dA"'{x) is the Euclidean volume element on the space 
of gauge potentials A and volQ^, is the volume of gauge transformations, 
viz., volume of 5'^7(iV)-valued functions on space. From (1) we see that 
SA = -D{5MM-'^), SA = D{M^-^SM'<) which gives 

dii{A) = (det DD) d^x{M, M^) (2) 

where dfj,{M,M^) is the volume for the complex matrices M,M^, which is 
associated with the metric dsj^ = 8/Tr(5MM-i Mt-i^Aft). This IS given 
by the highest order differential form dV as dji^M^ M^) = Y\^dV{M, M^). 
Writing out this form and using M = Up, H = = M'^M, 

dii{M,M'<) = l[dV{M,M^) vol{g,)^vol{g,) l[di^{H) (3) 

X X 

dp{H) = ea,...aAH-'dH)a,...{H-'dH),^ 
vol{g,) = l[dp{U) (4) 

X 

dp{U) is the standard Haar measure for SU{N). The volume element or 
the integration measure for the gauge-invariant configurations can now be 
written as 

dp{A) _ [dA^dA^] 



= {detD,D /^'^^/jf^ ^ {det DD)dp{H) (5) 

The problem is thus reduced to the calculation of the determinant of the 
two-dimensional operator DD. This is well known to be given in terms of 
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the Wess-Zumino-Witten (WZW) action as [6] 



det'dd 



n dimG 



(detDD) 



exp [2cA S{H)] 



(6) 



where caS"^ = it is equal to N for SU{N). The WZW action5(i/) 

for the hermitian matrix field H is given by [7] 



We can now write the inner product for states |1) and |2), represented by 
the wavefunctions and ^^2, as [8] 



The next step is the change of variables in the Hamiltonian. However, 
there is some further simplification we can do before taking up the Hamilto- 
nian. The wavefunctions, being gauge-invariant, are functionals of the matrix 
field H, but actually they can be taken as functionals of the current of the 
WZW model (7) given by J = {cA/'K)dzH H~^. Equation (8) shows that 
matrix elements of the theory are correlators of the hermitian WZW model 
of level number 2ca- The properties of the hermitian model of level number 
k + 2c A can be obtained by comparison with the SU (iV)-model defined by 
QkS{U)^ U(^£^ g SU{N). The hermitian analogue of the renormahzed level 
K — {k-\- ca) of the SU {N)-m.odie\ is —{k -\- ca)- Since the correlators involve 
only the rcnormalizcd level k, the correlators of the hermitian model (of level 
(A; + 2c^) ) can be obtained from the correlators of the SU (A^)-modcl (of level 
A; ) by the analytic continuation n —>■ —k. For the SU{N)f;-m.ode\ there are 
the so-called integrable representations whose highest weights are limited by 
k (spin < k/2 for SU{2), for example). Correlators involving the noninte- 
grable representations vanish. For the hermitian model the corresponding 
statement is that the correlators involving nonintegrablc representations are 
infinite. In our case. A; = 0, and we have only one integrable representation 
corresponding to the identity operator (and its current algebra descendents) . 
Therefore, for states of finite norm, it is sufficient to consider J [1, 2]. 




X 



Tr{H-^d^HH-^d^HH-^daH) 



(7) 




(8) 
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This means that we can transform the Hamiltonian Ti. = T + V to express 
it in terms of J and functional derivatives with respect to J. For functionals 
which can be expanded in powers of J, the kinetic term may worked out 
by the chain rule of differentiation, calculating the coefficients of derivatives 
with respect to the J's using proper regularization. We find 



T 
V 



1 

2? 



Ju 
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u 



+ 



TT 



mcA 



5J"-{u) 

dJa{x)dJa{x) 



{u) 5J\v) 



(9) 
(10) 



where m — b^ca/'^'t^ and 



ca 



TT^ (u 



7r(u 



(11) 



The first term in T shows that every power of J in the wavcfunction gives a 
value m to the energy, suggesting the existence of a mass gap. 

The next step is to solve the Schrodinger equation. We will consider the 
vacuum wavefunction, which is presumably the simplest to calculate. 

For the vacuum state, we take the ansatz — exp(P), where P is taken 
to be a series in powers of J. Substituting this into the Schrodinger equation, 
with the Hamiltonian given by (9, 10), P can be determined term by term. 
Upto three powers of the current we find 

P = / S'^(X1)X(X1,X2)5"(X2) + I r''(xi,X2,X3)r{xi)j''{x2)r{xs) + ■■■ 



K{xi,X2) 



[m 



V2 



(12) 



-I Xl,X2 



/ (xi, X2, X3) is given in reference [3]. The first term in (12) has the correct 
(perturbative) high momentum limit, viz., 

1 



exp' 



X) 



B^{y) + 0(3J) 



(13) 



The terms with higher number of J's can be shown to be small for the low 
momentum limit and for the high momentum limit. 
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Wc now use this result to calculate the expectation value of the Wilson 
loop operator which is given as 



W{C) = TrP §c^^dz+Adz) ^ rjYp g(T/cA)/c 



J 



(14) 



For the fundamental representation, its expectation value is given by 



where Ac is the area of the loop C. a is the string tension. This is a predic- 
tion of our analysis starting from first principles with no adjustable parame- 
ters. Notice that the dependence on and is in agreement with large- A'" 
expectations, with a depending only on the combination e'^N as iV — > oo. 
(The first correction to the large- A?" hmit is negative, viz., —{e^N)/2N'^\/^ 
which may be interesting in the context of large- analyses.) Formula (15) 
gives the values v^/e^ = 0.345,0.564,0.772,0.977 for A^ = 2,3,4,5. There 
are estimates for a based on Monte Carlo simulations of lattice gauge the- 
ory. The results for the gauge groups SU{2), SU{3), SU{4:) and SU{5) are 
Vo'/e^ = 0.335, 0.553, 0.758, 0.966 [9]. We see that our result agrees with 
the lattice result to within ~ 3%. 

One might wonder at this stage why the result is so good when wc have not 
included the 3 J- and higher terms in the the wavefunction. This is basically 
because the string tension is determined by large area loops and for these, it 
is the long distance part of the wavefunction which contributes significantly. 
In this limit, the 3J- and higher terms in (12) are small compared to the 
quadratic term. We expect their contribution to a to be small as well; this 
is currently under study. 

Another feature of the wavefunction is that P is nonlocal when expressed 
in terms of the magnetic field. This is essentially due to gauge invariance 
combined with our choice of = 0; it has recently been shown that a similar 
result holds for the Schwinger model [10]. 

Some observations on the magnetic mass of YM^^i at finite tempera- 
ture can now be made. The expression (9) shows that for a wavefunc- 
tion which is just J", we have the exact result T J" = m J". When 
the potential term is added, J" is no longer an exact eigenstate; we find 



{Wf{C)) 



constant exp [—(7^(7] 




(15) 
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(T + V) J"' = \Jm^ — J" + • ■ ■, showing that the mass value is cor- 
rected to the relativistic dispersion relation. J" may be considered as the 
gauge-invariant definition of the gluon. This result thus suggests a dynam- 
ical propagator mass m = e^c^/27r for the gluon. A different way to see 
this result is as follows. We can expand the matrix field J in powers of 
which parametrizes H, so that J ~ {cA/'K)dip,J'-a- This is like a perturbation 
expansion, but a rcsummcd or improved version of it, where we expand the 
WZW action in exp(2cA5(i7)) but not expand the exponential itself. The 
Hamiltonian can then be simplified as 



where (j)a{k) — y CAkk / {271111) (pa{k) , in momentum space. This result is 

obtained by expanding the currents and also absorbing the WZW-action 
part of the measure into the wavefunctions, i.e., the operator (16) acts on 
i|f _ gCAS(i?)^ ^]^g above equation shows that the propagating particles in 
the perturbative regime, where the power series expansion of the current is 
appropriate, have a mass m = c^ca/^tt. This value can therefore be iden- 
tified as the magnetic mass of the gluons as given by this nonperturbative 
analysis. 

For SU{2) our result is m ~ 0.32e^. Gauge-invariant resummations of 
perturbation theory have given the values 0.38e^ [11] and 0.28e^ [12]. Lattice 
estimates of this mass arc 0.31e^ to 0.40e^ (as a common factor mass for 
glueballs [13]) and 0.446^ to 0.466^ in different gauges [14] . An unambiguous 
lattice cvahiation of this would require a gauge-invariant definition a gluon 
operator and its correlator. Philipsen has recently given a definition of the 
"gluon" on the lattice. Preliminary numerical estimates then give the mass 



Let us now turn to excited states. Eventhough J is useful as a description 
of the gluon, it is not a physical state. This is because of an ambiguity in 
our parametrization (1). Notice that the matrices M and MV{z) both give 
the same A, A, where V{z) only depends on z and not z. Since we have the 
same potentials, physical results must be insensitive to this redundancy in 
the choice of M; in other words, physical wavefunctions must be invariant 
under M —>■ MV{z). J is not invariant; we need at least two J's to form an 




(16) 




as 0.376^ 



[15]. 



7 



invariant combination. An example is 



*2=/ fix,y) \dJa{x){H{x,y)H-\y,y))ab9My)] (17) 

Jx,y 



This is not quite an eigenstate of the Hamiltonian. Neglecting certain 0{J^)- 
terms, one can show that this is an approximate eigenstate of eigenvalue E 
if 



^Jm'^ - Vf + ym^^^vf + log(|a;-|/|7A) 



f{x,y) = E f{x,y) (18) 



where A is a scale parameter [2] . An a posteriori justification of this would 
require that the size of the bound state be not too large on the scale of 1/m. 
This does not seem to be the case, so, at least for now, all we can say is that 
the mass of ^2 > 2m; see however [2]. 



3 The Yang-Mills-Chern-Simons Theory 

We have extended our analysis to the Yang-Mills-Chern-Simons theory by 
adding a level k Chern-Simons term to the action, which gives a perturbative 
mass e^/An to the gluon [16]. The inner product now becomes 

(1|2) = J dii{H)e^^+^'^^ "^^^^ ^1^2 (19) 

The earlier conformal field theory argument shows that there are now new 
integrable operators. These lead to screening behaviour for the Wilson loop 
operator rather than confinement, as expected. The kinetic energy term 
becomes 

The gluon mass is now {k + 2cA)e^/47r reflecting dynamical mass generation 
as well by nonperturbative effects. The vacuum state and a number of excited 
states have also been constructed in this case [16]. 
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4 Towards the 3 + 1 dimensional theory 



The success of our analysis for YM2+1 depended largely on the use of the 
paramctrization (1). This facilitated the calculation of the inner product and 
the change of variables in the Hamiltonian. S first step in attempting 

a generalization to FM3+1, we can try to get a suitable paramctrization for 
a gauge potential A", i — 1,2,3, on R^. {Aq can be taken zero as before.) 
Recently, A. Yelnikov and I have found a paramctrization which may be 
appropriate [4]. For an U{N) gauge field, the matrix M in terms of which 
the potentials can be parametrized must take values in U{2N, C); this leads 
to more redundancy than in the 2 + 1 dimensional case and will require 
further constraints. Specifically, let be a basis of the Lie lagebra of U (N) 
(considered as NxN matrices), and let CTj denote the Pauli matrices. Then a 
basis for U{2N, C) is given by {1 1 ai®t"-, ai®it°-}. One can then 
show that for every gauge potential ^4?, one can construct a U {2N, C)-valued 
matrix M by solving the equation 

a.A = -a- dM M'^ (21) 

Conversely, if an arbitrary M E U {2N, C) is given, one can obtain a gauge 
potential from the above equation, provided M obeys the further con- 
straints 

Tr(tV • dM M-^ - h.c.) = 

di{M^aiM) = (22) 

These results hold for a neighbourhood of the fiat potential in the space of 
potentials. If we expand M as M 1 + if^ip"- + iakt"'Q% + • • •, it is easy to 
see that these reproduce known Abelian gauge field parametrizations. In an 
NxN splitting, we may write 

M=f^ 0) N (23) 



. [/. 

where can be obtained as an algebraic function of the gauge- invariant 
quantities H = MW and Wi = M^aiM. U can be removed by a gauge 
choice. The invariant measure of integration is then given by 

= d/.(i/, M^) exp(r + f ) 

5[a ■ dNN-^ + h.c] dlTxf'a ■ dNN-^ - h.c] (24) 
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dfi{H, W) is the product over the spatial points of the volume on U {2N, C)/U{N) 
and can be explicitly given [4] . Unlike the lower dimensional case, F can only 
be calculated in powers of F^; to the lowest order it is given by 



We are now trying to work out the Hamiltonian in these variables. The 
hope is that the use of F + F, upto the order calculated, will suffice to show 
that the wavefunctions should contain a term like exp(— / F'^/S/i) with a 
dimensional parameter /i (related to the A parameter). In this case, by 
comparison with YM2+1, one can conclude that the string tension should be 
y/a = /j,^ {N"^ — l)/87r. This is under investigation. 
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